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ABSTRACT 

We present simulations of initially stable isothermal clouds exposed to ionising 
radiation from a discrete external source, and identify the conditions that lead to 
radiatively driven implosion and star formation. We use the Smoothed Particle 
Hydrodynamics code SEREN and an HEALPix-based photo-ionisation algorithm 
to simulate the propagation of the ionising radiation and the resulting dynamical 
evolution of the cloud. We find that the incident ionising flux, ^Lyc critical 
parameter determining the cloud evolution. At moderate fluxes, a large fraction 
of the cloud mass is converted into stars. As the flux is increased, the fraction of 
the cloud mass that is converted into stars and the mean masses of the individual 
stars both decrease. Very high fluxes simply disperse the cloud. Newly-formed 
stars tend to be concentrated along the central axis of the cloud (i.e. the axis 
pointing in the direction of the incident flux). For given cloud parameters, the 
time, t^, at which star formation starts is proportional to '^^^ pattern 

of star formation found in the simulations is similar to that observed in bright- 
rimmed clouds. 

Subject headings: hydrodynamics — methods: numerical — stars: formation — H II 
regions 
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Introduction 



When an expanding Hll region overruns a pre-existing cloud, it drives an ionisation 



front and a shock wave into the cloud (Sandford et ah 1982 Bertoldi 1989 Lefloch & 



LazarefF 1994). As a consequence, the inner parts of the cloud are compressed, and may 
become gravitationally unstable, collapsing to form new stars. At the same time, the outer 
parts are "boiled off" by the ionisation front, and the rest of the neutral gas is stretched out 
into a cometary globule, with a dense head surrounding the newly-formed stars, a bright 
rim pointing towards the ionising star, and a cometary tail pointing away from the ionising 
star. 



There are many observations suggesting that star formation has been triggered inside 



cometary globules and bright-rimmed clouds (e.g. Lefloch & Lazareff 1995; Lefloch et al. 



1997 


Sugitani et al. 


1999 


2000 


2002; 


Walborn et al. 


2002 


Morgan et al. 


2008 


Smith et 



al. 2010). However, the extent to which star formation has actually been triggered by the 



shock front that preceeds a D-type ionisation front - as opposed to simply being revealed 



by the dispersal of residual gas - is not always clear (e.g. Dale et al. 2007 Indebetouw et 



al. 2007 Guarcello et al. 2010 Elmegreen 2011). There is some evidence that the young 



stellar objects deeply embedded inside a bright-rimmed cloud are the youngest ones, whilst 
those close to the rim of the cloud or distributed inside the expanding Hll region are 



older (e.g. Ikeda et al. 2008; Chauhan et al. 2009); this implies sequential triggering. In 



addition, Sugitani et al. (1999, 2000) report that the young stellar objects detected inside 
bright-rimmed clouds tend to lie close to the line joining the centre of mass of the cloud to 
the ionising star. 

Simulations of the interaction of ionising radiation with self-gravitating clouds have 



been presented by various authors. Kessel-Deynet & Burkert (2003) show that including 



self-gravity leads to greater compression, and may therefore trigger star formation. 
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Gritschneder et al. (2009) show that stable Bonnor-Ebert spheres can be driven into 



gravitational collapse, and that the final mass and age of the collapsed core depend on the 



incident fiux of ionising radiation. Miao et al. (2009) find that the morphological evolution 



of a cloud is very sensitive to the relative importance of self-gravity. Gritschneder et al. 



(2010) and Mackey & Lim (2010) have evaluated the roles played by turbulence and thermal 



instability in the formation of pillar-like structures. Henney et al. (2009) and Arthur et al. 



(2011) have explored the way in which magnetic fields infiuence the evolution of an Hll 



region, and conclude that they may inhibit fragmentation of the neutral gas swept up by 
the expanding Hll region; their simulations include both a proper treatment of the thermal 
microphysics and detailed predictions of diagnostic line emission, but not self-gravity, so 



they are focussing on different issues from us. As noted by Deharveng et al. (2005) "no 



model explains where star formation takes place (in the core or at its periphery) or when 
(during the maximum compression phase, or earlier)". 

In this paper we perform simulations of pre-existing clouds exposed to ionising radiation 
from a discrete source, and identify the circumstances that lead to triggered star formation. 
The paper is organized as follows. In Section |2] we define the initial conditions. In Section 
[3] we describe the numerical methods. In Section |4] we present and analyse the results. We 
summarize our main conclusions in Section |5l 



2. Initial conditions 

We define the initial conditions with reference to two coordinate systems, as illustrated 
schematically in Figure [1} 

The first co-ordinate system is centred on the ionising star. Distances from the ionising 
star are designated by the variable D. In particular, the position of the ionisation front in 



D IF. 



Fig. 1. — Configuration of a cloud of radius 



that is ionised by an external star 



(solid black dot). The cloud is located at distance -Dcloud fro^i the star. The dashed curve 
inside the cloud demarks the position of the ionisation front (I.F.) at the beginning of the 
simulation. 

the direction of the unit vector h is given by D^^ (fi) n, and the centre of the cloud is initially 
at position -Dcloud where k is a unit vector along the z-axis (see Figure [l]). The ionising 
star is characterised by the rate, A/'l^c, at which it emits hydrogen-ionising photons. 

The second co-ordinate system is centred on the initial centre of mass of the cloud. 
Distances from this point are designated by the variable R. The fiducial cloud is modeled 
as a Bonnor-Ebert sphere, i.e. an equilibrium isothermal sphere ( Bonnor||1956 Ebert 1957). 
It can therefore be characterised by its mass, M^lqud — ^ ^0 ' isothermal sound speed, 
Oq = 0.2kms~^ (corresponding to molecular gas with fraction by mass of hydrogen X = 0.7 
and temperature T = lOK), and its Bonnor-Ebert parameter, C.^^^^ = 4 (the dimensionless 
radius). By choosing < 6.451 we ensure that the cloud is initially stable, and therefore 
that its subsequent collapse is due to the incident ionising radiation. With this choice of 
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parameters, the fiducial cloud initially has central density and radius given by 

(^be) in-20 -3 /i\ 

^CENTRE = 4/^3^2 ^ °gCm3, (1) 

CLOUD 

D _ ^ ^CLOUD ^BE ^ nSnr f?! 

where /^(^be) ^'^^ dimensionless mass obtained by solving the Isothermal Equation 
( Chandrasekhar| 1939 Chandrasekhar & Wares 1949). 



The ionisation routine is designed to treat the D-type expansion of an Hll region. 
Therefore each simulation starts when the ionisation front switches from i?-type to D-type 



(Kahn 1954). We assume that at this juncture the matter has not had time to move. On 



Figure [T] the ionisation front is denoted by a dashed line. 

We have performed three suites of simulations. The first suite explores the consequences 
of changing Af^^y,^, and hence the flux of ionising photons, ^^yc ~ -^Lyc /^^-^cloud ' i^icident 
on the fiducial cloud. In the second suite, we explore the consequences of changing Mq^oud; 
keeping a^, and -Rcloud /-^cloud fixed, with a view to determining the maximum ionising 
flux that triggers star formation. In the third suite we explore the effect of numerical 
noise, by invoking different random - but smooth - initial distributions of SPH particles. 
In all simulations, the cloud is placed at distance -Dcloud — fro^i the ionising 

star, so that the divergence of the ionising radiation impinging on the cloud, measured by 
^cloud/^cloud = 0.1, is invariant. 

3. Numerical methods 
3.1. SPH Code 



We use the Smoothed Particle Hydrodynamics (SPH) code SEREN, which has been 
designed to investigate star and planet formation problems, and extensively tested on 
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a wide range of standard problems, as described in Rubber et al. (2011) and on the 
SEREN webpag^ In the present work we use the following SEREN options: standard 
SPH ( Monaghan] [1992 ) with a fixed (and exact) number of neighbours J^j^sEm — ^Oj 
standard artificial viscosity with a^y = 1, = 2 and no switches; an algorithm based on 
HEALPifl for following ionising radiation (see Section [sTsj below; and |Bisbas et al. | |2009[ ); 



and a barotropic equation of state for the neutral gas (see Section 3.4 below). The SPH 
equations of motion are solved with a second-order Leapfrog integrator, in conjunction 
with an hierarchical block time-stepping scheme. Gravitational forces are calculated using 



an octal-spatial tree (Barnes & Hut 1986), with monopole and quadrupole terms and a 
standard geometric opening-angle criterion. The tree is also used to generate neighbour 
lists. 



3.2. Initial equilibrium cloud 

The initial cloud configuration is created by first settling a cube of SPH particles of size 
2x2x2, using hydrostatic forces and periodic boundary conditions, but no self-gravity, to 
obtain a uniform-density glass-like distribution of particles. Next a unit-radius sphere is cut 
from this cube, the number of particles in the sphere, A/'gpjjgpj,, is counted, and the particles 
are given equal masses m^p^ = ^cloud /-^sphere- Finally the particle positions are stretched 
radially to reproduce a Bonnor-Ebert density profile; for a particle initially at radius r, we 
solve the equation 

= ^VI^be) (3) 



"^http: / /www. astro.group.shef.ac.uk/seren 
^http: / /healpix.jpl. nasa.gov/ 
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for ^, and then displace the particle radially according to 



-^CLOUD ^ 



(4) 



In all simulations we arrange that the mass of an SPH particle is m^p^^ ^5x10 ^ Mg . 
Given the barotropic equation of state that we invoke (Eqn. [6]), the minimum Jeans mass is 



^MiN ~ 0.005 Mg, and is therefore just resolved with ~2A/'i^ 



100 SPH particles (Bate 



& Burkert 


1997 


Whitworth 


1998 



3.3. Propagation of the ionising radiation 



The ionisation routine uses the HEALPix algorithm (Gorski et al. 2005) to create an 
hierarchy of rays emanating radially from the ionising star (cf. Abel & Wandelt 2002). 
In the immediate vicinity of the star there are only twelve rays, each at the centre of 
an approximately square element of solid angle = Att/12 1 steradian. Each ray is 
then split adaptively into 4 child-rays, and the splitting is repeated recursively so that the 
separation between neighbouring rays is always less than one half of the local smoothing 
length (i.e. always less than one eighth of the diameter of the local SPH particles). Thus, 
after i splittings each ray is at the centre of an approximately square element of solid angle 
= steradian, and is separated from its neighbours by an angle 9^ ~ 2~^ radian. In the 
present simulations there are up to eleven levels of splitting, so the rays on the finest level 
are ~ 1.5arcmin apart. At each step we randomise the orientation of the twelve rays on 
level (and hence also the orientation of rays on higher levels), in order to avoid numerical 



artefacts (cf. Krumholz et al. 2007). 



Along a ray with direction given by the unit vector fi, the ionisation front is located at 
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position Djp(n) n, where 



p\Dh)D^dD = ^^^^ 



(5) 



D=0 



Here p{Dh) is the mass-density at position Dn, m = rUp/X is the mean mass associated 
with each hydrogen nucleus, rUp is the proton mass, X = 0.7 is the fraction by mass of 
hydrogen, and is the recombination coefficient into excited stages only (i.e. we invoke 
the On-The-Spot Approximation; Osterbrock 1974). The integral in Eqn. ([s]) is evaluated 
using a second-order integration scheme, with an integration step that is one quarter of the 
local smoothing length (i.e. one sixteenth of the diameter of the local SPH particles). 



Further details of the method are given in Bisbas et al. (2009). 
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Fig. 2. — False-colour maps of the initial gas temperature on slices through the middle of 
the fiducial cloud, for three different ionising fluxes. From left to right, ^^^(. = 10^cm~^s~^ 
(Simulation 1), ^^^^ = 3 x 10^°cm"^s"^ (Simulation 4), and ^^yc = lO^^cm^^s"^ (Simula- 
tion 7). The colour bar gives log^g(T/K), and the axes are in pc. The snapshots are taken 
after the flrst iteration of the code; the time on each map denotes to the moment when the 
ionisation front switches from R to D type. 
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3.4. Thermodynamics 

Neutral gas. The temperature, T^^, of the neutral gas is calculated using a barotropic 
equation of state, 

P ^ 



TM=To{l+[-^] \. (6) 

V P CRIT / 

where = 10 K, p^mt — 10^^'^ g cm~'^ and 7 = 5/3 is the ratio of specific heats. This 
equation of state is presumed to mimic the gross behaviour of protostellar gas: at low 
densities, p ^ Pcrit; S^^ tends to be optically thin to its own cooling radiation and 
therefore approximately isothermal at T ~ 10 K; at high densities, p ^ Pcrit; ^^e gas tends 
to be optically thick to its own cooling radiation and therefore heats up adiabatically. Since 
the gas dynamics is only followed to densities ~ lO^^^gcm"^ (thereafter sink particles are 
introduced; see below), the temperature is always T < 200 K. Therefore we can assume 
that 7 ~ 5/3, even for H2, because the rotational degrees of freedom of H2 are only weakly 
excited. 

Ionised gas. The temperature, Tj, of the ionised gas is assumed to be Tj = 10*^ K, except 
in the immediate vicinity of the ionisation front, where the temperature changes smoothly 
from Tj to T^^ over a region with a width of two smoothing lengths (i.e. half the diameter of 



an SPH particle; see Bisbas et al. 2009). Figure |2j shows initial temperature cross-sections 
through the fiducial cloud for three representative ionising fluxes, at the time when the 
ionisation front switches from R to D type. 



3.5. Sink particles 



We follow Bate et al. (1995) in replacing the densest regions with sink particles. 
Specifically, if the density of SPH particle i exceeds Pgi^K = lO^^^gcm"^, then SPH 



particle i and its neighbours are excised from the simulation and replaced by a sink particle 
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having the same mass, centre of mass and momentum as the sum of the SPH particles it 
replaces. Any SPH particle which subsequently passes within -Rgjj^j^ — 4AU of a sink, and 
is gravitationally bound to it, is accreted by the sink, i.e. it is excised from the simulation 
and its mass and momentum are assimilated by the sink. The presumption is that the 
material in a sink will inevitably collapse to form a single star or a close multiple system. 
We therefore hereafter refer to sinks as stars. However, we do not take account of radiative 
feedback from these newly-formed stars. 

The value of Pgi^K is chosen so that any protostellar condensation is already well into 
its Kelvin- Helmholtz contraction phase (and therefore well relaxed), but not yet into the 
second collapse phase, when it is replaced with a sink. With Pgjf^j^ = lO^^^gcm"^, the 
temperature inside a protostellar condensation has increased from ~ lOK to ~ 200 K, by 
the time it becomes a sink. The sink radius is not specified, but is calculated on-the-fly 
to ensure that a protostellar condensation that becomes a sink is resolved; this gives 

-"-SINK — V'-'^"min/ ^" rsiNK/ — *±X^W. 



3.6. Magnetic fields 



The SEREN code is not able to handle magnetic fields, and so our results must 
be interpreted as exploring what can happen when the magnetic field is dynamically 
unimportant. It will be important to explore this aspect of the problem, but even if 
a uniform magnetic field were assumed, it would necessarily introduce two further free 
parameters (field strength and direction). Previous studies of irradiated clouds embedded 
in a uniform magnetic field suggest that a sufficiently strong magnetic field may seriously 



affect the cloud evolution (WiUiams 2007 Henney et al. 2009; Mackey & Lim 2010). 



Depending on the initial field orientation with respect to the direction of ionising photons, 
the cloud is either fiattened (perpendicular orientation) or the radiative implosion may be 
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prevented (parallel orientation). However, these findings were not confirmed by Arthur et 



al. (2011) who studied formation and evolution of globules at the border of the Hll region 



expanding into the turbulent, magnetised molecular cloud, and found that the magnetic 
effects are weak. Further studies are needed to resolve this discrepancy. A turbulent field 
and/or non-ideal effects would complicate matters still further. 



4. Results 

4.1. The fiducial cloud exposed to different ionising fluxes 

In the first suite of simulations, we treat the fiducial cloud (with mass M^^^^j^ = 5 , 
isothermal sound speed = 0.2kms~^ and Bonnor-Ebert parameter = 4, hence radius 
-^cLouD ~ 0.3 pc; we place it at distance D^^^^^ = lOR^^^^^ = 3pc from the ionising star; 
and we simply vary the ionising output of the star, M^^^, and hence the flux of ionising 
photons incident on the cloud, $Lyc- Table[l]gives the ID, A/'l^^ and ^^yc ^"^^ each simulation 
in this suite. 

Table 1: The ID, output of ionising photons from the central star, M-^^^-, and ionising flux 
incident on the cloud, for the simulations involving the fiducial cloud. 



Simulation ID 


1 


2 


3 


4 


5 


6 


7 8 




1048 


3 X 10^« 


1049 


3 X 10^9 


1050 


3 X 10^° 


10^1 3 X 10^1 


/cm"^s"^ 


10^ 


3 X 10^ 


10^0 


3 X 10^° 


10" 


3 X 10" 


10^2 3 X 10^2 



4- 1.1. Overview 



In all simulations, the ionising radiation erodes the cloud from one side, and 
the ionised gas flows away into the surrounding space. There is a range of ionising 
fluxes, that trigger star formation in the cloud; for the fiducial cloud this range is 
10^cm~^s~^ < $^^^(2 ^ 3 X 10^^cm~^s~^. If the ionising flux is very low, the ionisation 
front is R-type, and the cloud is just slowly eroded. If the ionising flux is > 10^cm~^s~^, 
the ionisation front becomes D-type, and is preceded by a shock front, which propagates 
into the cloud compressing it. There is then a competition between erosion of the cloud 
by the ionisation front and collapse of the cloud triggered by the shock front. Provided 
the ionising flux is not too large, ^^^^ ^ 3 x 10^^ cm~^ s~^, at least some of the cloud ends 
up in newly-formed stars. However, if the ionising flux is very large, the cloud is rapidly 
dispersed, with no star formation. All simulations are followed until all the cloud gas is 
either in newly-formed stars, or has been ionised and dispersed. The time at which the first 
star forms is denoted by . We limit further discussion to the cases in which star formation 
occurs, namely those with IDs 1 to 6 in Table [Tj 



4- 1.2. Phenomenology of star formation for low ionising fluxes 

If the incident ionising flux is at the low end of the range for triggered star formation, 
the shock front driven into the exposed side of the cloud is weak, and it sweeps up matter 
rather slowly. Consequently, long before the shock front reaches the centre of the cloud, 
the pancake-shaped layer of swept-up gas becomes sufficiently massive to contract laterally 
due to its own self-gravity. As the material in the layer converges on the axis of symmetry, 
the increase in pressure deflects the inward flow towards the centre of the cloud, creating 
a dense filament down the axis of symmetry. It is in this filament that star formation first 
occurs, usually at a point well behind the ionisation front and towards the centre of the 
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cloud. For these relatively low ionising fluxes, the ionisation front on the exposed side of the 
cloud develops a U-shape structure. Figure [3] is a time sequence of column-density images 
for Simulation 1, illustrating these morphological features, and terminating at the moment 
the flrst star forms. 



4.1.3. Phenomenology of star formation for high ionising fluxes 

If the incident ionising flux is at the high end of the range for triggered star formation, 
the shock front driven into the exposed side of the cloud is strong, and it sweeps up matter 
rather rapidly. As a result the swept-up layer becomes sufficiently massive to fragment, 
before it has had time to undergo much lateral contraction under its own self-gravity, and 
therefore before there has been time to form an axial filament. Consequently, star formation 
first occurs near the centre of the swept-up layer, and close to the ionisation front. For 
these relatively high ionising fiuxes, the ionisation front on the exposed side of the cloud 
develops a V-shape structure. Figure |4] is a time sequence of column-density images for 
Simulation 6, illustrating these morphological features, and terminating at the moment the 
first star forms. 



4.1.4. Formation time of first star 

Figure [5] shows column-density images of all the simulations involving the fiducial cloud 
and resulting in star formation, at the moment star formation occurs, t^. This sequence 
demonstrates how, as the ionising fiux is increased, the length of the central filament 
shortens, the location of star formation moves closer to the ionisation front, and the shape 
of the ionisation front morphs from a broad U-shape to a narrower V-shape. As noted by 
Gritschneder et al. (2009), star formation occurs earlier for higher ionising fiuxes. Figure |6] 



is a plot of against $Lyc, for the fiducial cloud. The dashed line shows the numerical 
results, and the solid line shows the best power-law fit to these results, 



4.1.5. Cloud morphology at 

The morphological evolution of the cloud and the location of star formation can be 
quantified coarsely by defining two metrics at (see lefthand panel of Figure [?]). 

The first metric is the width, of the cloud at t^, i.e. the maximum extent of 
the neutral gas perpendicular to the direction of the ionising star. The variation of 
with log^g($Lyc) is shown in the middle panel of Figure [?[ On average, decreases with 
increasing 3>LyC' froi^ ^* > 0.3 pc at the low fluxes corresponding to broad U-shaped 
ionisation fronts, to < 0.2 pc at the high fluxes corresponding to narrow V-shaped 
ionisation fronts. 

The second metric is the distance, from the ionisation front to the first star, at 
t^. The variation of d^ with log^jj(<l'Lyc) shown in the righthand panel of Figure [7j d^ 
decreases monotonically with increasing $LyC' i^om d^ > 0.08 pc at low $LyC' ^* < 0.02 pc 
at high $ 



4.1.6. Star formation after t^ 

After t^, further star formation is usually confined to the region between the first star 
and the ionisation front. Once the shock front has passed, any neutral gas that has not yet 
been incorporated into stars is overpressured, and tends to re-expand, as noted previously 



by Sandford et al. (1982); Bertoldi (1989); Lefioch & Lazareff (1994); Gritschneder et al. 
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(2009 


); 


Miao et al. 


(2009 



For low ionising fluxes, the expanding neutral material forms a single massive pillar 
protruding into the Hll region, with a dense head pointing towards the ionising star; this is 
where most of the additional stars are formed. Because of their location, they are quickly 
overrun by the ionisation front. This terminates accretion, and therefore they tend to have 
low masses. The external appearance of these stars as they emerge into the Hll region is 



very reminiscent of the Evaporating Gaseous Globules (EGGs) seen in M16 (Hester et al 



1996) and elsewhere. 



Figure |8] illustrates the remnants of the single pillar formed in Simulation 2, at 
t = 0.68 Myr. The three stars marked A have already been overrun by the ionisation front, 
and so their accretion has been terminated. The close binary system marked B is an EGG, 
just emerging from the head of its nascent pillar; its accretion rate is rapidly declining. The 
star marked C is still embedded in the pillar, and growing by accretion; at the end, it is the 
most massive star. 

For higher ionising fluxes, a fragmented bundle of smaller pillars is formed. The 
individual pillars do not contain sufficient mass to spawn many additional stars, and they 
are quickly dispersed by the ionisation front. 



4- 1-7. Star formation efficiency 



Figure [9] shows that the total mass of stars, M^^^, formed from the fiducial cloud 



is a monotonically decreasing function of the ionising fiux, $Lyc- Figure 10 shows that a 
range of stellar masses is formed, from 0.01 to 1.44 M^, with the upper envelope on the 
individual masses being a decreasing function of ^^yc • 

When the ionising fiux is low, the evolution of the cloud is relatively slow. Consequently 
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there is more time for stars to form, and - once they have formed - to accrete, because 
they remain embedded longer. As a resuh the total mass converted into stars is large 
(~ 0.8M|2j^Qyjj in Simulation 1), and at the end there are some relatively high-mass stars. 

When the ionising flux is high, the evolution of the cloud is relatively rapid, and a large 
fraction of the mass is eroded by the ionisation front before it can become dense enough to 
collapse and form stars. Moreover, stars that do form enjoy only a brief embedded phase 
during which they can grow by accretion. The total mass converted into stars is therefore 
small (~ 0.06M|_,j^Qyj5 in Simulation 6), and the individual stellar masses are low. 

There are too few simulations here to infer more detailed constraints on the statistics 
of stars formed by radiatively driven implosion, but extended protostellar discs, with radii 



> 100 AU, appear to be common, and there are several multiple systems. Figure 11 shows 
a face-on projection of a typical protostellar disc, taken from Simulation 3. The star is 
represented by the white dot with diameter 10 AU, at the centre of the disc — although the 
corresponding sink is actually half this size, with diameter 5 AU. At this juncture, the star 
is in a wide binary system with separation ~ 1000 AU. This binary subsequently hardens, 
due to dissipation between the circumstellar discs attending its individual components, and 
at the same time it moves towards a third star and forms an hierarchical triple. 



4.2. Non-fiducial clouds exposed to different ionising fluxes 

In the second suite of simulations, we treat clouds with different masses from the 
fiducial cloud, namely M^^^^^^ = 2.5 M^, 10 and 15 M^, but the same isothermal 
sound speed, = 0.2kms~^, the same Bonnor-Ebert parameter, = 4, and the same 
flux divergence -Rcloud /-^cloud — O-^' hence radii -Rcloud — 0.15pc, 0.6 pc and 0.9 pc and 
distances -Dcloud — 1-5 pc, 6pc and 9pc, respectively. Here we are only concerned with 
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locating the separatrix between ionising fluxes that do trigger star formation, and those 
that do not. Therefore, for each cloud mass we only report the results for two ionising fluxes 
(one either side of the separatrix), and we terminate the simulations either at t^, when the 
first star forms, or when only 20% of the cloud mass remains. Table [2] gives the ID, M^^^^^, 
-^cloud; -^cloud; -^Lyc "^Lyc simulations in this suite (Simulations 9 to 14). 



Figure [12] shows, for each of four cloud masses, pairs of ionising fluxes separated by 0.5 
dex, for which the lower flux (marked with a cross) triggers star formation and the upper 
one (marked with an open triangle) does not, but simply disperses the cloud. This plot 
shows that it requires a higher flux to disperse a less massive cloud without also triggering 
star formation. This is because the emission measure through a Bonnor-Ebert sphere 
with fixed isothermal sound speed and fixed Bonnor-Ebert parameter, is proportional to 
"^cLouD' Consequently, a less massive cloud is denser, and only a very high ionising flux 
can penetrate right through it before the ionisation front switches to D-type and starts 
compressing it. 



In support of this contention, we also plot on Figure 12 the critical ionising flux for 
which the ionisation front has just reached the central, approximately uniform-density 
part of the cloud (^ < 1) when it switches to D-type. This is obtained by evaluating the 
emission-measure integral 

m = f e-'^^^'U^', (8 



where ip{$,) = \'D.{p/ p^^^^^^^^) is the isothermal function ( Chandrasekhar fc Wares|[l949 ). The 
critical ionising flux is then 

^ ^ «B^^V^(eBE) [((^be) - C(l)] 

-> 6 X lO^^^cm'^s-^ ( ^m"""^ ) ' 
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where the final expression (Eqn. 10) has been evaluated for the fiducial parameters 

= 0.2kms~^ and = 4. For ionising fiuxes in excess of 
dispersed, without any star formation being triggered. 



the cloud is simply 



4.3. Influence of numerical noise 



In the third suite of simulations, we treat three additional realisations of the simulation 
involving the fiducial cloud and an ionising fiux ^^^^^ = 3 x 10^° cm~^ s^^. These simulations 



are identical to Simulation 4 (presented in Section 4.1), except that in each case the 
distribution of SPH particles representing the cloud has been rotated through three random 
angles. Consequently, the particle noise - which infiuences the details of when and where 
star formation occurs, and therefore what the final stellar masses are - is different. 



Table |3] summarises the global metrics for Simulations 4, 4a, 4b and 4c; and Figure 13 
shows surface-density images for the same simulations at t^, = 0.06 Myr (top row) and at 
t = 0.5 Myr (bottom row). All the global metrics (t^, d^,, Wi, and M^^^) are broadly similar, 
as is the morphology of the cloud at and the location of star formation. However, as 
would be expected for a chaotic non-linear process, the details of the final state (masses 
and locations of individual stars at t = 0.5 Myr) show a large spread, and this should be 
interpreted as a measure of the intrinsic variance of the process. 



5. Conclusions 

We present the results of simulations of star formation triggered by radiatively 
driven implosion, i.e. stable clouds which are driven into collapse and fragmentation by 
being overrun by an Hll region. The clouds are all modelled as stable Bonnor Ebert 
spheres (i.e. equilibrium isothermal spheres) with isothermal sound speed = 0.2kms~ 
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Table 2: The ID, cloud mass, M^^oud; cloud radius, -Rcloud; distance from ionising star to 
cloud, -Dcloud; output of ionising photons, A/'l^c; cind incident ionising flux, $LyC' 



non-flducial clouds simulated in order to generate Figure 12 
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Table 3: The ID; time at which star formation starts, t^; lateral width of the cloud, w^, 
when star formation starts; distance behind the ionsation front, d^, at which star formation 
flrst occurs; and total mass of stars formed, M^^^, for Simulation 4 and the three additional 
simulations (4a,4b,4c) performed with different realisations of the fiducial cloud exposed to 
ionising flux ^j^^^^ = 3 x 10^° cm~^ s~^. 

Simulation ID 4 4a 4b 4c 
0.06 0.055 0.06 0.06 
0.2 0.23 0.23 0.22 
0.029 0.034 0.024 0.027 
1.84 1.93 1.61 2.12 



4/pc 
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and Bonnor-Ebert parameter = 4. The ionising star is always placed at distance 
-^cLouD — froHi tlie Initial centre of the cloud, but its ionising output, Af^^y,^, and 

hence the ionising flux incident on the cloud, $LyC' varied through over three orders of 
magnitude. 

There is a range of ionising fluxes that trigger star formation. In this range, the time at 
which the flrst star forms is given approximately by = 0.19 Myr ($j^y^/10^ cm~^ s~^) 

For ionising fluxes at the low end of this range, the evolution of the cloud is quite 
slow. When star formation starts at t^, the ionisation front on the exposed side of the 
cloud has a broad U-shape, and star formation is concentrated in a dense fllament along 
the axis of symmetry of the cloud, some distance behind the ionisation front. This flnding 
is in agreement with observations of star formation in bright rimmed clouds, as reported 
by Sugitani et al. (1999, [2000) ). Although compression of the cloud is relatively weak, star 



formation is quite efficient (a large fraction of the cloud mass is converted into stars) and 
the individual stars have relatively high mass, because there is a long period during which 
star formation can occur; also, once formed, the stars remain embedded in the neutral gas 
for a long time and can therefore continue to grow by accretion. During the later stages of 
cloud dispersal, the remnants of the cloud appear as a pillar protruding into the Hll region; 
the newly-formed stars emerge from this pillar, and into the Hll region, as Evaporating 



Gaseous Globules (EGGs; see Hester et al. 1996). 



For ionising fluxes at the high end of this range, the evolution of the cloud is more 
rapid. When star formation starts at t^, the ionisation front on the exposed side of the 
cloud has a narrow V-shape, and star formation is concentrated at the tip of the V, just 
behind the ionisation front. Although compression of the cloud is relatively strong, star 
formation is quite inefficient (a small fraction of the cloud mass is converted into stars) and 
the individual stars have relatively low mass, because there is very limited time available 
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for star formation before the cloud is dispersed by being ionised; also, once formed, the 
stars do not remain embedded in the neutral gas for long, and so they cannot grow much by 
accretion. During the later stages of cloud dispersal, the cloud breaks up into a bundle of 
small pillars which do not have enough mass to spawn further stars. Once the stars emerge 
into the Hll region, their growth ceases. 
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density plots were made using the SPLASH visualization code (Price 2007). 
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Fig. 3. — Time sequence of surface-density images showing the evolution of the fiducial cloud 
when it is exposed to a relatively low ionising flux (Simulation 1, ^^^^^^ = 10^ cm~^ s~^). Note 
the overall broad U-shape of the ionisation front bounding the exposed side of the cloud, and 
the dense filament forming down the centre of the cloud. t/Myr is given in the top righthand 
corner of each frame. The colour bar gives log^g(S/gcm~^), where S is the surface-density, 
and the axes are in pc. 
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Fig. 4. — Time sequence of surface-density images showing the evolution of the fiducial cloud 
when it is exposed to a relatively high ionising flux (Simulation 6, ^^^^ = 3 x 10^^ cm~^ s^^). 
Note the overall narrow V-shape of the ionisation front bounding the exposed side of the 
cloud, and the sharp density peak at the apex of the cloud pointing towards the ionising 
star. t/Myr is given in the top righthand corner of each frame. The colour bar gives 
log (S/gcm~^), where S is the surface-density, and the axes are in pc. 
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Fig. 5. — Surface-density images of the fiducial cloud at the moment the first star forms, t^, 
for different ionising fluxes, as given in Table [T] top row, left to right. Simulations 1, 2 and 
3; bottom row, left to right. Simulations 4, 5 and 6. t^/Myr is given in the top right corner 
of each frame. The location of the first star is marked with a dark blue dot. The colour bar 
gives log^|j(S/gcm~^), where S is the surface-density, and the axes are in pc. 
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Fig. 6. — Logarithmic plot of the time at which the first star forms, against the ionising 
flux, ^. The solid line is the best fit power-law given by Eqn. (It)). 



Fig. 7. — Left panel: Snapshot of Simulation 1 at t^, illustrating schematically the definitions 
of and d^. Middle panel: as a function of log^g($Lyc) Simulations 1 to 6. Right 
panel: d^ as a function of log^^{^^^^) for Simulations 1 to 6. 
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Fig. 8. — Surface-density image of Simulation 2 at t = 0.68 Myr. 
log^jj(S/gcm~^), where S is the surface-density, and the axes are in 
for discussion. 



The colour bar gives 
pc. See Section 4.1.6 
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Fig. 9. — Log-linear plot of the ionising flux, ^^^^^ against the total mass of stars formed, 
Mj,Q^, from the fiducial cloud in Simulations 1 to 6. 
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Fig. 10. — Log-log plot of the ionising flux, ^^^^^ against the masses, M^, of the individual 
stars formed from the flducial cloud in Simulations 1 to 6. For Simulation 4 we have also 
plotted - at slightly displaced positions - the masses of the stars formed in the three addi- 
tional simulations performed to evaluate the influence of numerical noise (Simulations 4a, 4b 



and 4c; see Section 4.3); because of the intrinsic variance in the outcomes of these different 
realisations, this produces a small bulge in the upper envelope on the individual masses. 
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Fig. 11. — Surface-density image of a protostellar disc, taken from Simulation 3 at time 

for commentary. The colour bar gives log (S/gcm~^), 



t = 0.17 Myr. See Section 



4.1.7 



where E is the surface-density, and the axes are in pc. 
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Fig. 12. — Log-linear plot of ionising flux, ^^y,^ against cloud mass, M^^^^^ locating the 
separatrix between configurations that produce stars (SF) and those that do not (no-SF). 
All other critical parameters are invariant, viz. isothermal sound speed in the cloud, = 
0.2kms~^, Bonnor-Ebert parameter, ^g^, = 4, and flux divergence, -Rcloud /-^cloud — O-l- 
The crosses mark configurations that result in star formation, and the open triangles mark 
configurations that do not. The solid line shows the locus predicted by Eqn. 
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Fig. 13. — Surface-density images of Simulations 4, 4a, 4b and 4c. The first row is at 
= 0.06 Myr when star formation starts. The cloud morphology is broadly similar in all 
four simulations, and star formation is located at approximately the same distance ci^ behind 
the ionisation front. The second row is at t = 0.5 Myr, and the distribution of newly-formed 
stars and residual neutral gas is more chaotic, but statistically similar. The colour bar gives 
log (S/gcm^^), where S is the surface-density, and the axes are in pc. 
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